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Abstract. Ultrasound modulated optical tomography, also called acousto-optics tomogra- 
phy, is a hybrid imaging modality that aims to combine the high contrast of optical waves 
with the high resolution of ultrasound. We follow the model of the influence of ultrasound 
modulation on the light intensity measurements developed in [14]. We present sufficient con- 
ditions ensuring that the absorption and diffusion coefficients modeling light propagation can 
locally be uniquely and stably reconstructed from the corresponding available information. 
We present an iterative procedure to solve such a problem based on the analysis of linear 
elliptic systems of redundant partial differential equations. 

1. Introduction 

This paper concerns the reconstruction of the diffusion and absorption coefficients 7 and cr 
for a model of hght propagation in tissues given by 

-V • 7VU + au = in X 
u = f on oX, 



where u{x) is the hght intensity of diffuse photons propagating in a bounded open domain 
X C for n > 2 with Dirichlet boundary condition given by / on (the sufficiently smooth) 
dX, the boundary of X. More precise models of boundary conditions can be used with no 
consequence on the results presented in this paper and so we present our result with Dirichlet 
conditions to simplify notation. 

The reconstruction of (7, a) from optical boundary measurements is known to be a severely 
ill-posed problem in many settings of light propagation and results in reconstructions with 
[very poor resolution [6, 9]. Ultrasound Modulated Optical Tomography (UMOT), also called 
'Acousto-Optics Tomography (AOT), aims to combine the high contrast of optical coefficients 
.observed in many diseases with the high resolution of ultrasound. We refer to [7, 21, 30] for 
• • 'additional information on this hybrid modality. In this paper, we follow the model of ultra- 
. ^ -sound modulation of light intensities developed in [14]. In this model, the difference of light 
^ 'measurements with and without ultrasound modulation provides, to a first approximation, 
^ 'internal functionals of the unknown parameters of the form 

" " " (2) H{x) = 7(x)|Vn|^(a;) + r]a{x)u^{x), x £ X, 

where is a known constant in the model. The internal functional H{x) is parameterized 
by the boundary condition / on dX. The objective of UMOT is to reconstruct the unknown 
parameters (7, a) from knowledge of a minimum of functionals Hj for 1 < j < J corresponding 
to well-chosen boundary conditions {fj}i<j<j- 

Many studies are devoted to the problem with a = above. This problem, which finds 
applications in ultrasound modulated electrical impedance tomography (UMEIT), also called 
acousto-electric tomography or impedance-acoustic tomography, is now well understood and 
we refer the reader to [3, 11, 12, 13, 16, 18, 22, 24, 25, 26] for a detailed account of such theories. 
Note that [24, 25] analyze the reconstruction of tensor-valued (anisotropic) coefficients 7, which 
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we do not consider in this paper for UMOT. The reconstruction of two coefficients with more 
measurements than proposed above is considered in [4]. 

Pioneered by the work in [16], several papers [11, 12, 24, 25, 26] provide exphcit recon- 
struction procedures for 7 (when a = 0) in the setting of a highly redundant measurements 
corresponding to a large value of J. Such procedures do not extend to the reconstruction of 
(7, a) from knowledge of Hj in (2) for 1 < j < J. Rather, we follow a standard method to 
solve the above nonlinear problem consisting of analyzing its linearization and introducing a 
standard fixed point iterative procedure provided the linearized operator is injective. Such 
a strategy was followed in the setting cr = in, e.g., [8, 22, 23]. The papers [8, 23] provide 
general strategies to solve similar problems. In this paper, we follow the method developed 
in [8] that recasts the UMOT inverse problem as an elliptic redundant system of partial dif- 
ferential equations for which we can apply the classical elliptic regularity results developed in 
[1, 17, 29]. 

More precisely, we recast the UMOT inverse problem as a system of partial differential 
equations in section 2 and present some sufficient conditions so that the system be elliptic. 
Ellipticity provides an inversion procedure up to the possible existence of a finite dimensional 
kernel. This is not sufficient, and conditions for injectivity are given here in section 3. Once 
the linearized UMOT inverse problem is injective, standard fixed point iteration procedures 
recalled in section 4 may be applied to solve the nonlinear problem locally. The conditions 
for ellipticity and injectivity are shown to be satisfied for appropriate choices of the boundary 
conditions {fj} in section 5. The practically and pedagogically interesting case of constant 
backgrounds is presented in section 6. 

The UMOT inverse problem is an example of a large class of hybrid inverse problems (also 
called coupled-physics or multi-wave inverse problems) that aim to combine one modality with 
high contrast with another modality with high resolution. For an incomplete list of books and 
reviews on this active field of research, we refer the reader to, e.g., [2, 5, 10, 27, 28]. 



2. System and ellipticity 

We consider the equation 

—V • 7Vn,- -|- cTiio =0 in X 

(3) ' 3 J 

Uj = fj on dX, 

for 1 < j < J with J > 2 and availability of functionals of the form 

(4) 7|Vnj|2 _^ wlujl"^ = Hj{x) in X, I < j < J, 

for some fixed, known, constant ^ t] £ M. The equations (3)- (4) may be seen as a redundant 
2J X (2 + J) system of nonlinear partial differential equations for the dependent variables 
(7, C7, {uj}i<j<j)- As indicated above, in the setting cr = 0, explicit reconstruction procedures 
for (7,{uj}) have been obtained in [11, 12, 16, 24, 25, 26] when J is sufficiently large. The 
extension of such methods to the reconstruction of (7, a) directly from knowledge of {-Hj}i<j<j 
is not known at present. Moreover, as shown in [8] in the setting a = 0, 7 can uniquely and 
stably be reconstructed from {Hj}i<j<j for a smaller value of J than what is necessary in the 
available explicit reconstruction procedures. 

The method followed in [8], which shares many similarities with that described in [23] 
(the main difference being that [23] considers the inversion of systems of pseudo-differential 
operators while [8] considers the inversion of larger, but differential, systems of operators), 
first consists of linearizing the nonlinear problem (3)-(4). Replacing 7—7-7-1- ^7, a ^ a + 6a, 
and Uj — )• Uj + 5uj as well as Hj — ?■ Hj -\- 6Hj and considering the terms in (4)- (3) that are 
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linear in "5", we calculate that 

(5) |Viijp(57 + r]\uj\'^6a + 2^\7uj ■ \76uj + 2r]aoUj5uj = SHj, 

(6) L6uj - V • 6-fVuj + Uj6a = 0, L := -V • 7V + a. 
Moreover, we obtain that 5uj = on dX. Let us introduce the notation 



Fi = Vui, 



iVujI' 



and the operator 
(7) 



Then, applying L to (5), Mj to (6), and introducing the commutator [L,Mj] = LMj — MjL, 
we find that 

(8) (^|i^jf + MjV ■ 6-fFj'^ + (^Lr]\uj\^6a - Mjujda^ + [L, Mj]6uj = LSHj, 

(9) L6uj - V • 6-f\/uj + Uj6a = 0. 

This is a linear system of equations for v = ((^7, 6a, {5^ij}i<j<j), which may be recast as 

(10) Ajv = {Pj + Qj)v = Sj, 

with Pj a 2 J X (2 + J) matrix of second-order operators given by 

/ ; ; ; ; 

-7|FjpA + 27Fj : V® V 77|njfA ... [L,Mj] ... 

... -7A ... 

V ^ ; ; : : y 



(11) 



pj 



and Qj a 2J x (2 + J) matrix of at most first-order differential operators, whose explicit 
expression we do not reproduce. The source 5j is a 2J x 1 matrix with odd entries 2j — 1 
given by L6Hj and vanishing even entries. Here aOa for an n— vector a is the rank-one matrix 
of components aiaj; and for two matrices A and B, we denote A : B = Tr{A*B) the usual 
inner product with A* , the Hermitian conjugate to A. 

The lower-order term Qj has a complicated structure, which we hope not to analyze in detail. 
We thus look for conditions (on the boundary conditions fj and on the number of measurements 
J ) which guarantee that Pj is an elliptic operator. Let pj{x,£,) be the (principal) symbol of 
Pj. The operator Pj is said to be elliptic when for each G S"""^, the 2 J x (2 -|- J) matrix 
Pji^: has full rank 2 + J. The latter matrix is given by 



(12) pj{x,0 



\ 



7(x)|F,|2(x)|ep-27(x)(F,(x)-0' -ri\u,\Hx)\e 




7(x)iep 



V 



/ 



with Pj{x, ^) a quadratic form in ^ whose explicit expression does not influence ellipticity since 
7(a;)|i^P is uniformly bounded from below by a constant times which equals 1 for ^ € S^~^. 
Therefore, pj{x,Q has full rank for a given ^ G if and only if the following J x 2 matrix 
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has full rank (i.e. rank equal to 2) 



(13) 



\ 



We assume, for instance by imposing that fj is uniformly bounded from below by a positive 
constant, that \uj\ is also uniformly bounded from below by a positive constant. As a conse- 
quence, we observe that the above matrix is full rank if and only if the following matrix is full 
rank 



(14) 



\ : ■■ J 

with the dependence in x dropped to simplify notation. We define the quadratic (in ^) forms 

(15) pj{x, = lei' - 2ej{x) ® 9j{x) : ^ ^3 e = icp - mi^) ■ 0'- 

Were we to be in the situation where the absorption coefficient a vanishes in the above equa- 
tions, then pj is full rank (i.e. of rank 1) provided that {pj{x,S^) =0, 1 < i < >/} implies 
that ^ = 0. In other words, the light cones with direction 6j intersect only at the point ^ = 0; 
see [8] for an analysis of this problem in the linearized setting. 

Here, in the general a case, we obtain that the matrix pj{x, ^) is of rank two for all ^ G S"~^ 
if and only if the quadratic forms 

(16) Pjk{x,0 = \dj\'^{x)pk{x,(,) - \dk\'^{x)pj{x,S), 
are such that 

(17) {Pjk{x,i) = ^, l<j<k<j} implies that ^ = 0. 

Let us assume that the ellipticity condition (17) holds. This provides uniform ellipticity 
inside the open domain X. We now need to find boundary conditions that preserve the 
elliptic structure of Pj. Such conditions are called the Lopatinskii conditions, and under 
the assumption (17), it is not difficult to show that they are satisfied for Dirichlet boundary 
conditions for v; see [8]. 

Thus, using the theories of elliptic systems described in [1, 29] (see also [8]), we obtain 
results that we collect in the following theorem: 

Theorem 2.1. Consider the system (10) augmented with boundary conditions v = on dX. 
We assume that Pj is elliptic in the sense that (17) holds. Then the system admits a left 
parametrix Bj such that 

(18) V = BjSj + Tjv, 

where Tj is a compact operator. More precisely, Tj is a pseudo- differential operator of order —1 
(of class L^Q as defined in, e.g., [19]J. The above expression may be extended by a parametrix 
that solves the above system up to a smoothing operator of arbitrary order. That is, we have 

m— 1 

(19) v=Y,rjBjSj + Tfv 

for any m where TJ^ is a pseudo- differential operator of order —m, mapping functions (or 
distributions) in H^{X) to functions in H'^^'^'{X). Moreover, there exist constants C and C2 
such that 

(20) \\{5-i,5a)\\Hs(x-M.'') + ¥uj\\Hs+i(^x-^-') < C'll^^j IIh^cx^^) + C2\\v\\L^(^X;^^+Jy 
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This estimate is optimal since \\SHj\\H''{X;R-^) bounded by a constant times the above right- 
hand side. 

Proof. Elliptic theory in [1, 29] provides a left parametrix Bj such that 
(21) Bj{Pj + Qj)v = {I- Tj)v = BjSj, 

where Tj is a compact operator, and more precisely a pseudo-differential operator of order 
— 1 (of class L^l as defined in, e.g., [19] when the coefficients are assumed to be sufficiently 
smooth). When the coefficients in the differential operators are sufficiently smooth, the per- 
turbation Tj may be replaced by an arbitrarily smoothing operator. That is, for any m we 
have 



(/ - T] 



m—l 



T'jBjSj, 

q=0 



where TJ" a pseudo-differential operator of order — m, bounded from H^{X) to H'^^'^{X). 
Moreover, standard elliptic regularity results [1, 29] provide the following stability estimates 



(22) 



\'^\\h=+'2{X;R2+J) < C\\Sj\\h=(^X;R'^J) ^ C!\\6Hj\\fjs+2(^X;RJ) + C2\\v\\l2(^X;R2+J)- 



These estimates are not optimal as the source term in (9) vanishes. From standard elliptic 
regularity for the latter equation, we obtain that Suj is one derivative smoother than ^7 and 
hence we obtain (20). That the latter estimate is optimal is clear from (5). □ 

Remark 2.2. The estimate (20) can be obtained directly from (8)-(9) by introducing an 
elliptic system in the sense of Doughs and Nirenberg [17]. Assigning the weights S2j-i = 
to (8) and S2j = — 1 to (9), and the weights ti = t2 = 2 for the columns for {6j,6a) and the 
weights tj = 3 for 3 < j < 2 + J for the columns for {6uj}, we find that the leading term in 
(8)-(9) is 

/ : : : : :\ 



(23) 



Pj 



-7|Fj|2A + 2-fFj (g)Fj :V (S)V r/|MjfA 









-7A 



V • : : : : y 

This operator is elliptic if and only if Pj is elliptic and standard weighted elliptic regularity 
estimates [1, 29] then directly provide (20); we refer to [8] for additional details, which we do 
not reproduce here. 

Remark 2.3. The same estimates may in fact be established at a lower cost in terms of 
the number of redundant boundary conditions. Indeed, consider the problem (5)-(6) directly. 
Then we can also write this system in the sense of Douglis and Nirenberg assigning the weights 
S2j-i = to (5) and S2j = — 1 to (6), as well as the weights ti = t2 = for the columns for 
((57, 5a) and the weights tj = 1 for 3 < j < 2 -|- J for the columns for {6uj}. For such weights, 
we would find that the leading term is of the form 

/ : : : : :\ 



(24) 



Vj 



-Fj-V 



r]\uj\^ 







2-/Fj • V 
-7A 



y : : : : :/ 

It is not difficult to verify that the symbol of this operator (replacing V by i^ and A by — 1^|^) 
is injective if and only if Pj above is elliptic. Moreover, we verify as in the case u = in 
[8] that under the same conditions of ellipticity, then the Lopatinskii conditions are satisfied 
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provided that 6uj = 0; with no conditions necessary for {5a, 5j) on dX. Elhptic regularity 
then again provides a result of the form (20) with the last term replaced by C2||(5uj||j;^2(j5f.]Rj). 
In some sense, this latter result is more favorable and shows that the application of L to (5) 
above requires that we introduce additional boundary conditions for (SajSj) on dX. Since 
the injectivity results of the following section can only be established for systems of the form 
(8)-(9), we decided to directly present the latter in detail. Nonetheless, the following more 
precise version of Theorem 2.1 holds: 

Theorem 2.4. Consider the system (5)-(6) with Dirichlet conditions 5uj = on dX . Then 
under the same ellipticity condition (17) as in Theorem 2.1, the existence of parametrices for 

V solution of (5)-(6) as described in (18) and (19) still holds. Moreover, we have the (optimal) 
stability estimate 

(25) ||('^7,^Cr)llH»(X;IR2) + \\SUj\\Hs+i^X;RJ) < C\\6Hj\\H^(X;RJ) + ^^2 || || ^2 (^.jj,/) . 

3. Generic injectivity 

The presence of the constant C2 in the preceding estimates comes from the fact that the op- 
erator Pj + Qj augmented with Dirichlet conditions need not be injective. This is reminiscent 
of the possible lack of invertibility of —A + V for a given potential V. Elliptic theory provides 
a parametrix A^^ (inversion of the Laplace operator with Dirichlet boundary conditions) so 
that —A + y is replaced by I — A^^V, with A^^y a compact operator as Tj above. However, 
we cannot guarantee that 1 is not an eigenvalue of Tj in general unless we have a precise 
understanding of the lower-order term Qj. 

We present here a general methodology to ensure generic injectivity of an operator of the 
form Pj+Qj at the cost of imposing additional boundary conditions. Again, consider the scalar 
problem Av := (—A-\-V)v = with prescribed Dirichlet conditions on dX. Let us replace the 
above system by A*Av = 0, which is a fourth-order system, with its own Dirichlet conditions, 
which in this case would correspond to prescribing v and d^v on dX. By integrations by parts, 
we obtain that A*Av = along with v = d^v = on dX implies that Av = with the same 
boundary conditions. When V is real-analytic, then the Holmgren theory for Cauchy problems 
can be invoked to show that v = 0, yielding injectivity for the non homogeneous problem. 

If A is injective, then so is A + 6V for 6V sufficiently small since A~^SV is of spectral radius 
strictly less than 1 for 6V sufficiently small (in an appropriate topology and for an open set 
that depends on V). As a consequence, we obtain that there is a dense open set of potentials 

V (in that same appropriate topology) such that A + V is injective. This property (valid on 
a dense open set) is called generic injectivity. 

When V is not analytic, then a considerably more delicate Unique Continuation Property 
(UCP) for A shows again that V = 0. We do not consider UCP here and refer to [8] for the 
analysis of the case a = 0. 

Coming back to the UMOT problem, (10) says that Ajv = Sj, which we can modify as 

(26) A*jAjv = A*jSj in X, v = d^v = g on dX. 

Since we prescribe Dirichlet data in (1) to be the same as that of the background solutions, 
f = on the boundary. However, the data for d^Uj is then determined by these solutions 
of (1), and therefore d^duj is nonzero in general. Nonetheless, we do assume this is known 
boundary data and represent it as g in (26). 

Theorem 3.1. [Holmgren] Let us assume that all coefficients (7,cr, {uj}) are real-analytic so 
that A*jAj is an operator with real-analytic coefficients. Assume moreover that Pj is elliptic 
in the sense that (17) holds. Then AjA j and Aj augmented with the boundary conditions 
v = d^v = on dX are injective operators. Moreover, the problem (26) has a unique solution 
satisfying 

(27) II^^IIh»(X,R2+J) < C \ \\5Hj\\Hs(^x,'&J) + IIS'll//s-3/2(0X;R2+./) j • 
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Proof. Injectivity is a direct consequence of Holmgren's theory for systems of operators as 
presented, e.g., in [20] (see also [8, Thm 3.3]). The function g G H''-^/'^{dX,R'^+-^) can be 
lifted to the unique M^+'^) function (j) such that 

(28) AA0 = in X; cj) = 0, d^(t> = g on dX. 
By ellipticity, we can apply (22) to u — </> to obtain 

(29) \\v\\H''iX,R'2+J) < C (\\6Hj\\Hs(^x,RJ) + \\4>\\h''(X;R'^+J)) + C'z || L2(x.ir2+ J) 

and injectivity precisely means that C2 = up to the choice of a possibly larger constant C. 
We also know that 

(30) II</'IIh={v,rj+2) < C'lbll//=~3/2(av,RJ+2) 

from which the result follows. □ 

With this result, we obtain a generic injectivity statement for the linearized UMOT problem. 

Theorem 3.2. [Generic Injectivity] There is a dense open set of coefficients (7,0") (in any 
topology of sufficiently smooth coefficients) such that for J > 5 (J > 3 in dimension n = 2) 
and for an open set of boundary conditions {fj}i<j<j , we have that A*jAj and Aj augmented 
with the boundary conditions v = dyV = on dX are infective operators, and that the problem 
(26) has a unique solution. 

Proof. Let (7,0") be real-analytic. Then the results of Theorems 5.1 and 5.2 below show 
the existence of an open set (in the C^{dX) topology, say) of boundary conditions {/,} such 
that Pj is elliptic. By density, we choose the boundary conditions {fj} to be real-analytic. As 
a consequence, the solutions Uj are real-analytic and Pj is an elliptic system of second-order 
operators. 

We may then apply Theorem 3.1 to obtain that AjAj and Aj are injective. Now, injectivity 
of A*jAj and Aj still holds for {'y,d) sufficiently close to (7,0") (for instance in the C^{X) 
topology) and for fj close to fj (for instance in the C^{dX) topology). □ 

The above construction thus shows that for appropriate choices of boundary conditions, the 
operator Aj is injective for an open, dense, set of coefficients (7,0"). Note that an open dense 
set may in fact be very far from accounting for all possible coefficients (7,0"). But this result 
still indicates that injectivity of Aj and A*jAj does hold for a large class of coefficients. 

Another injectivity result shows that when Pj is elliptic, then AjAj in (26) is injective 
on sufficiently small open domains $7 9 0. For x £ X, we say that x G eil. for < e < 1 if 
£-^x £ n. The size of the domain depends on the coefficients in Pj, but this results gives yet 
another indication that the UMOT problem is well posed in several configurations of interest. 

Theorem 3.3. [Injectivity for Small Domains] Let Q be a subset of X with G eQ for all 
< £ < Eq. Assume that the highest order part of the system Pj given by (11) is elliptic in 
the sense that (17) holds. Then there exists ei > such that if e < ei, the linearized system 
for V = (57, 6a, 6ui , . . . Suj)^ 

A*jAjv = A*jSj, on eQ. 
dv 

V = — = g on (9(ef]) 
dv 

has a unique solution v satisfying 

(31) \\v\\Hl(tn;RJ+'^) < Ci\\5H\\h2^^q.^J-) + Ce||5'||^l/2(g(^n).iRj+2) 

where Ci is independent of e. 

Proof. Recall that Aj = Pj + Qj where all terms of Qj are at most first order. We let 
A j(0) = Pj{0) + Qj{0) be the operator obtained by freezing all coefficients (7, a, {uj}, {Vuj}) 
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at X = 0. Then Aj = Aj{0) + [Aj - Aj{0)) and we set the remainder T = [Aj - Aj{0)) . Let 
us hft g to the unique function (j) defined by (28). Then we have that 

(32) \\4>\\H^en) < C'e|bllHl/2(a(en)) 

where may depend on e. Let us subtract off (j) so that w = v — (p solves 

A*jAjw = A*j{Sj + Aj^), on eQ 

w = -5^ = on d(en). 
ov 

Note that Pj{Q) also satisfies the ellipticity condition (17) by assumption, and is homoge- 
neous of degree 2. Since the coefficients of -Pj(O) are all constant, it is easy to see by using 
Plancherel's theorem that 

(33) / \\Pmwf>a\\w\\H2^^^y 

where a does not depend on e. Suppose 

A*j{^)Aj{{))w = B*f 

for w G HQ^eVi) where S is a second order differential operator. Then B is bounded indepen- 
dently of e from to by ||-B||. Integration by parts yields 

f Aj{f))w Aj{{))w = [ fBw 

JeQ Jen 

from which we can compute 

(34) [ Pj{0)w ■ PjiO)w = -2 [ QjiO)w ■ PjiO)w - [ Qj{0)w ■ QjiO)w + [ fBw 
Jn Jn Jn Jn 

(35) < C\\w\\Hi{en)\\w\\H2{en) + \\B\\\\f\\L^(en)\\'w\\H^{en) 

where we use the full and norms here. Given the Poincare inequality for the fixed 
domain Q and w E Hq{Q), 

we can deduce from scaling the following Poincare inequalities for the small domain eO, 

(36) \Mm(eQ) < D\\w\\h2^^^^, 

and 

(37) \Mm{en) < ^ellw^llnKen) 

for all w G HQ{eQ,). Applying (36) and (33) on the left hand side of (35) yields 

(^IMn^ien) ^ '^^ {\\M\m(en)\\w\\Hl{en) + \\B\\\\f\\L-^{en)\\w\\Hl{en)^ 
Now we use (37) and assume that e is small enough so that 

a - CD^e > 0. 

Bringing the first term on the right hand side over to the left and dividing by ||f ||j:/2(^f^) gives 
us 

CD 

(38) \\M\Him< ^_^j^, m\\f\\L^.n). 

This shows that {A*j{Q)Aj{Q))-^B* (with zero boundary conditions) exists and is bounded as 
an operator from Lp' to iJg, independently of small domain size e. Indeed we have 

CD 

(39) \\{AmAmr'B*\\ < —-^pn. 
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Now, for smooth enough coefficients, the remainder T = (^Aj — Aj{0)) will be small in norm 
from i/^(eO) to L^(eil) due to the smallness of the domain. We now want to show that 

(40) A*jAj = {Aj{0) + Ty{Aj{0) +T)w = 0, w = d,w = 0, 

implies that tf; = 0, in other words that (^j(O) +T)*(Aj(0) +T) is injective on Hq. The above 
implies that 

(41) w = {A}mj{0))-'{{Aj{0)+TnAj{0)+T)-A*jiO)Aj{0))w 

(42) = {A*j{0)Aj{0))-^{A*j{0)T + T*Aj{0) +T*T)w. 

Now, for e sufficiently small and the coefficients in the equations sufficiently smooth, T is small 
in norm from to L^. Hence applying (39) to each of the three terms in (42) implies that 
w = 0. Furthermore, the linearized inverse with the data subtracted off, w = v — (j) satisfies 

(^j(O) + T)*(^j(0) + T)w = {Aj{Q) + TfiSj + Aj(t>), 

from which the smallness of T, the fact that the components of S are either zero or L6Hj, and 
the estimate (39) will yield 

(43) \\M\Hl{en;RJ+2) < C (^\\SH\\ ^2 (^^q.^J-) + \\4>\\H2{en;RJ+2)) 

which, using (32), implies (31). □ 

We recall that in all the settings where AjA j can be proven to be injective, then we have 
the ellipticity estimate (22) with C2 = for w = v — (j) where (p solves (28). In other words, 

(44) 11(^7, 5cr) 11^3 (X;R2) + \\Uj\\Hsl^X;RJ) < C (^\\SHj\\Hs(^X;RJ) + ll5'll//^-3/2(aX;R-^+2)) ' 

This shows that the reconstruction of (^7, 6a) from knowledge of 6Hj and g is obtained with 
no loss of derivatives when Pj is elliptic and A*jAj is injective. 

4. Local reconstruction of the nonlinear problem 

The elliptic stability estimate (44) obtained in the preceding section when Pj is elliptic and 
A*jAj is injective holds for coefficients {'y,a,Uj) e X' = H'{X) x H'{X) x H'+'^{X;R'^) for 
s > ^. The availability of such an estimate allows us to solve the nonlinear problem (3)-(4) 
locally after a few algebraic manipulations. 

Let us recast the nonlinear problem (3)- (4) as 

(45) T2j^iiv) = 0, F2j{v) = Hj, l<j<J. 

The construction of the linear operator Aj in the preceding section is based on the commuta- 
tion relation (8). We need to construct an inverse problem for v with Aj as its Frechet deriva- 
tive. We also need to transform the above over-determined system into a determined system of 
equations if we want to invert it numerically. To do so, we consider a point vq = (7, a, uj) G Af* 
and assume that the linear differential operator Aj constructed in the preceding section with 
coefficients given by vq is such that AjAj augmented with boundary conditions as in (26) 
is injective. Let L and Mj be the operators in (6)- (7) obtained with the coefficients vq. We 
recast the inverse problem (45) as 

(46) J2j_i (t') = 0, J^2j (v) := LF2j {v) - MjF2j-i{v) = LHj := Kj, l<j< J. 
The above problem is recast using the more compact notation 

(47) Fiv) = K, 

with /C := LT-L the 2J-vector of sources with J of them non-vanishing. We therefore view the 
sources m a space identified with y^""^ = i?*~2(X; M"'). We also identify % with the sources 
{Hj}i<j<J in y' = H%X-RJ). 
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The above nonhnear problem is still over-determined and may not admit solutions unless /C 
satisfies compatibility conditions. We thus apply the adjoint operator Aj as we did in earlier 
sections and denote the solution of 

(48) A*j A jw = A*jS in X, w = dyW = goudX, 

as w = QqS + RoQ, with Qq mapping continuously 3^*"^ to X'^ and Rq mapping continuously 
i?^-i(aX;M2+J) to X'. 

Two nonlinear procedures may then be set up depending on the construction of vq. We can 
first assume that vq is such that J-'(fo) = /Co = {{LHj^o}). This may be constructed by using 
a first guess for (7,0") and then solving for the uj. We then observe that 

J^'{vo){v - Vo) = T{v) - F{vq) - {F{v) - F{vq) - F'{vq){v - Vq)) = K, - + Q {v - Vq), 

where G{v — vq) is quadratic m. v — vq and bounded in y^""^ by C\\v — folios • Applying A*j, 
we obtain 

A*jF'{vq){v - vo) = A*j{lC- /Co) + A*jg{v - vo), 

since J^'{vo) = Aj by construction of F in (46)-(47). Let us define g = dyV and = d^VQ on 
dX. Then we find that 

(49) v-vo = QoA*j{JC - /Co) + Ro{9 - 9o) + {A*jAj)^'A*jg{v - vo) = fo + :S{v - vo), 

with fo = QoA*j{K: - /Co) + Ro{g - go) and 5(0) = {A*jAj)^^A*jg{cl)). Here, {A*jAj)^^ is the 
inversion of the operator A*jAj with vanishing boundary conditions (of the form w = and 
dyV = on dX). 

For /o sufficiently small in , which means for % — Hq sufficiently small in 3^^* and g — go 
sufficiently small in H^~^ {dX; M^+'^), we then easily verify that (j)^ fo +Z{4>) is a contraction 
on a sufficiently small ball in X^ . This shows that the above equation (49) admits a unique 
solution that can be computed by the converging algorithm v^~^^ — vo = fo +^{v'^ — vo)- 

Alternatively, we may construct vo as the solution to 

A*jT{vo) = AjlCo in X, dyVo = g on dX, 

where we also prescribe vo on dX. The main difference with respect to the preceding setting is 
that here vo is constructed as the solution of a fourth-order system with its Neumann condition 
may be chosen so that d^vo = d^v, which we need to assume is known to obtain an invertible 
operator A*jAj. In such a setting, we obtain 

(50) v-vo = QoA*j{lC - /Co) + {A:jAj)-^^A*jg{v - vo) = fo+2{v - vo), 

with fo = QoA*j{n - Ho) provided that A*jAj remains invertible. For H — Tio sufficiently 
small in 3^**, we again obtain that the solution to the above equation is unique and computed 
by the converging algorithm v^^^ — Vo = fo + 2{v^ — vo)- 

In the latter setting, assume that v and v are two solutions of (47) with right-hand side 
given by /C = LT-L and /C = LT-L, respectively. Then we find that 

(51) II7 - 7||h»(x) + \W - ^Wn^ix) < CWH - 'H||^.(X;rj)- 

In other words, the iterative algorithm (50) converges to a solution that satisfies the optimal 
(elliptic) stability estimate (51). This result is summarized as: 

Theorem 4.1. Let vq = (co, 70, {lijo}) ^ be a point such that Pj is elliptic and A*jAj is 
injective. Let Ho he given by (47). Let v = {a,^,{uj}) G X and v = (a,7,{uj}) S X and let 
% and % he given in (47) hy solving the prohlems (3) -(4). Then ifH. and % are sufficiently 
close to T-Lq in y, then the error estimate (51) holds. 
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The above theorem, which is the main contribution of this paper, shows that the UMOT 
problem is well-posed in the sense that smah errors in the available internal functionals (in a 
norm H^{X) with s sufficiently large) translate into small errors in the reconstruction of (7, a) 
in the same norm. Moreover, the iterative algorithm based on (50) converges to the unique 
solution of the UMOT problem. In the setting of (49), we also obtain a converging algorithm 
with a similar estimate as (51) with an additional term of the form \\g — g\\ ^^_z{dX;'R.'^^'^) 
on the right-hand-side. As in the analysis of many hybrid inverse problems, our proof for such 
results requires that the coefficients one wishes to reconstruct be sufficiently regular [10]. 

5. Conditions for ellipticity. 

In this section, we collect several sufficient conditions that ensure the ellipticity of the 
operator Pj. We show that in the general case, there exist sufficient background solutions so 
that the linearized operator Pj is elliptic. Recall that 



and the system is elliptic if 

(52) m''-2ie,-if)d] = i\9,\^-2i9,-if)dj Vz,j = l...J 

has no solutions | = ||| e S"~^. 

In dimension n > 3, the following result says that locally four solutions are sufficient to 
guarantee ellipticity. This number climbs to five if we want to guarantee ellipticity on a domain 
of arbitrary size. The proof is based on the construction of complex geometric optics (CGO) 
solutions, see (53) below. The use of CGO solutions to solve such problems was introduced in 
[15]; see [10] for a review of the method. 

Theorem 5.1. Let X C M" a bounded open domain. Then there is an open set of boundary 
conditions {fj} for 1 < j < 5 such that the operator Pj is elliptic. 

Proof. Consider solutions of (3) of the form 

(53) ^^(^) = _LeP-(l + V>^(:c)), 

where p € is a complex vector such that p ■ p = 0. Note that the real and imaginary parts 
of Up above are then also solutions of the equation (3) with real-valued coefficients. 

It is shown in [15] that ipp is of order J/o]"^, e.g., in the C^{X) norm provided that the 
coefficients (7, a) are sufficiently smooth. As a consequence, all properties that can be shown 
when Ipp is set to 0, such as the absence of solutions to (52), extend to solutions of the form 
(53). We will therefore find five real- valued solutions of the form (53) such that (52) admits no 
solution. By continuity, any choice of fj sufficiently close (in any sufficiently strong topology, 
say C^{dX)) to the traces of said solution on dX, will generate solutions such that (52) admits 
no solution. 

It thus remains to prove that we can find five (locally four) solutions of the form J}?(e''''^) 
or $5(6*^'^') such that (52) admits no solution (we easily verify that for p sufficiently large, the 
presence of the term ^ does not modify our conclusions, see [10]). 

Let M be a large real number and p = \k\{ei + 262). We denote by ui and U2 the real and 
imaginary parts of UMp, that is 

Ml = e^'^l'^l^'i cos M\k\x2, U2 = e^^l'^'^'i sin M\k\x2 

We find 

cosMlA;l2;2 \ q - i sinMlA;l2;2 \ ^ _ cosMlA;l3;2 ^ _ s\D.M\k\x2 
— sinM]A;]x2 / \ cos 
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We also define U3 = ui + n2, which is also clearly a harmonic solution. Finally, we define M4 
and U5 exactly as ui and U2 above with M replaced by 1. For 1 < j < 3 and 4 < m < 5, we 
have from (52) 

{l-2{9,.0')dl = {l-2{e^-0^)dl 

Since d1 + = \k\~'^, at least one of them does not vanish (this is why we need five solutions 
whereas four solutions suffice locally with either d^ or d^ bounded away from 0). For M 
sufficiently large, we thus find that (l — 2{9j ■ ,^)^) is as small as we please for 1 < j < 3. But 
it is proved in [8] that this is not possible for any ^ G In other words, the latter results 

states that when a = 0, then three solutions such as ui, U2, and U3 above are sufficient to 
provide ellipticity. This concludes the proof of the theorem. □ 

In dimension n = 2, the above proof combined with the results in [8] may be used to 
show that four solutions (J = 4) guarantee ellipticity. However, in fact three well chosen 
solutions (via well-chosen boundary conditions fj) guarantee ellipticity as the following result 
demonstrates. 

Theorem 5.2. In dimension n = 2, the results of Theorem 5.1 hold for J = 3. 

Proof. Consider the same solutions Uj for 1 < j < 3 as in the preceding theorem with 
M = 1. Knowledge of Hj for 1 < j < 3 is equivalent to knowledge of 

-f^l2 = 7VU1 • VU2 + 7]aUiU2. 

This is clear by polarization and is a consequence of the fact that our internal functionals are 
quadratic in the solutions Uj. 

We then calculate that the operator Pj is elliptic, or equivalently (52) admits no solution, 
when the system 

- 2{0^. if )di = (1021'- 2(62 ■ if )di 
-201 ■ i02 ■ id\ = (1^1 12 - 2(^1 • if)did2 

-201 ■ i02 • idl = {\02\^ - 2(02 • if)did2 

has no solutions for unit ^ G S^. Note that since 0i and 02 are unit and orthogonal, 

01 • C = cos a, 02 • ^ = sin a 

for a the angle between 0i and ^. Let us assume that all the above equations hold. Plugging 
this into the above system means that 

(54) (1 - 2cos2a)(i^ = {1 - 2sm'^ a)dl 

(55) — 2cosasinadf = (1 — 2 cos^ a)(ii(i2 

(56) — 2cosQsinQ(i| = (1 — 2 sin^ a)(ii(i2- 
Equation (54) is recast as 

/ ■ 2 2 Nj2 / ■ 2 2 Nj2 

(sm Q — cos a)d2 = — (sm a — cos a)di 

which implies that 

sin a — cos a = 0, 

since di,d2 are nonnegative and not both zero. Now, either (55) or (56) implies that 

sin a cos a = 

which is a contradiction since we must have sin a = it cos a. So the symbol of the system is 
elliptic. The rest of the proof is the same as that of Theorem 5.1. □ 
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6. Constant Background. 

For the important case when the background 7, a are constant, we can actually back- 
substitute for 5uj in (5), (6) and still obtain a system of differential equations if one chooses 
the boundary data correctly. We this approach, we still get full ellipticity, injectivity and 
stability for the linearized problem, which here becomes an explicit 2x2 fourth order elliptic 
system of PDEs with constant coefficients. Let us take special background real CGO-like 
solutions i = 1 ... J of the form 



-X-Vi 



where the vj are vectors of unit length. With these (exact) solutions we have 

Oi = -7= — r = Vi 

\\/Ui\ 

and 

di 



|Vnj| V <^ 

Again we take data of the form Hi (4). Then all determinants of the system (13) vanish when 

for alH, J = 1 . . . J . Since we always have here that = and \6i\^ = 1, this is equivalent to 
(57) {i.v,f = {i.v,f 

for all So, we have ellipticity if we take enough Vj so that there is no unit length ^ for 
which this holds. When n = 2, it suffices to take the three background solutions with 

vi = ei,V2 = e2,v^ = (^/2/2,^/2/2). 

When the dimension n = 3, we can take the 4 directions 



Vi 



l<i<3, Vi = (l/\/3,l/^,l/\/3). 



For any dimension we can find such directions to guarantee ellipticity if we allow ourselves n+1 
background solutions (which is optimal in dimensions n = 2 and n = 3 to obtain ellipticity as 
demonstrated earlier). We can, for example take 

Vj=ej, 1 < j < n, Vn+l = {0,Wi, . . . ,Wn-l), 

where w € M"~^ is any unit vector for which no possible sums of the components of the form 

±Wi lb 7i;2 ± . . . lb Wn-1 = 1- 

Now, we know from using Holmgren's theorem in Section 3, the linearized problem here is 
invertible using boundary data described above. However, in this section we process the data 
slightly differently. Note that the ellipticity conditions don't change when we eliminate the 
unknowns 5uj from the system (5), (6). From (5), 

(58) SHi{d-f, 6a) = d-flVui]"^ + 2yV6ui • Vui + 7]6a{uif + 27]a5uiUi. 



Using 



we have that 



Vui = , l—UiVi, 
7 



(59) - — - = 5y—Ui + 2j. f^V5ui ■ Vi + rjdaui + 2ria6ui. 
Ui 7 V 7 
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Now, to eliminate Suj we apply L to both sides: 

(60) L — - = L{5^)—Ui S^Ui — aV5j ■ ViUi + 27. [—L(V6ui 

Ui 1 1 V 7 



r]L{5a)ui — ar]6aui — ^rjS/da ■ ViUi + 2r]aL{5ui), 



and we can calculate 



L{6ui) = ■ Vi. —Ui + S^—Ui - 5aui 



7 7 



which yields 

L(V6ui ■ Vi) = VL{5ui) ■ Vi 



We obtain 



V(V(57 • Vi) ■ Vi^^Ui + 2VSj ■ Vi^Ui + dj^^^Ui - SaJ'^Ui - V5a ■ ViUi. 



where Cj and are the second order constant coefficient differential operators 



and 



Define the data as 



a = - A + 2dl^ + (3 + 27?) J -a,, + 2(1 + 7?) ^ 



Bi = -77^A - (2 + Vi)J ^dn, - 2(1 + 7]). 
a \ a 



(61) Si = L5Hi = ( —L— ) 5Hi 



and define the operator 
(62) Aj = {C, Bi 



with as many rows as we have i's. This is a constant coefficient elliptic PDE system, but 
due to the possible presence of eigenvalues, Dirichlet conditions are not obviously enough to 
guarantee injectivity. To resolve this, we again apply the technique described in Section 3. 
That is, we make the problem square by applying the adjoint. In this case, however, the result 
is a 2 X 2 fourth order system of differential equations 

More precisely, we have 

A*jAj = 




\ i i / 

Since A*jAj is elliptic in the sense of Douglis-Nirenberg, the system will admit a parametrix if 
we apply the Dirichlet conditions w = Q and dyW = on dX . Again by using Holmgren, we 
have that the operator A*jAj is actually injective with these conditions. Indeed, let us assume 
that AyAjw = in X and w = and duW = on dX . Then integrating by parts, we obtain 
that Ajw = in X. Since Aj is elliptic, this implies that d'^w = so that all second-order 
derivatives of w vanish. Upon differentiating the equation for w, we obtain that all derivatives 
of order less than or equal to 3 vanish on dX (as well as higher-order derivatives in fact) . From 
Holmgren we deduce that w = in the vicinity of dX and by induction in the whole domain 
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X. So, the leading order operators are invertible stably for constant (7, u), and by continuity 
for background coefficients which are near to constant. We have proven the following theorem. 

Theorem 6.1. Suppose the background coefficients (7,0") are constant and the domain X C 
M" . Then if one applies data fi corresponding to n + 1 background solutions 



-X-Vi 

for the {vi} any set of unit vectors for which 

(63) (i.v,f = ii.v,f y{i,j) 

has no solutions ^ G S"""*^, then the linearized system for w = {6j, 6aY 

A*jAjw = A*jS in X 
dw 

w = — — = on dX 
ou 

for Aj and S given by (62), (61) is elliptic and has a unique solution w satisfying 

(64) WMlH^iX-X^) < C\\dH\\Hs(^X;RJ)- 

Remark. Note that in the above analysis we did not explicitly cover when a = 0. This case 
is actually simpler than the above, and illuminates very clearly why more than two solutions 
are necessary for ellipticity. We can of course choose the background solutions uq = 1 and 
Uj = Xj for 1 < j < n. We then obtain that 

dHo = r]6a. 

This is easy enough to invert for 6a, which we eliminate from the system of equations and 
then focus on 6^. After elimination, we find that 

AdHi5^ = {A - 2dl)5^ 

which is a hyperbolic operator (in fact it is precisely the wave operator where Xi acts as time). 
So, if we were only to use one Uj, (in addition to no), it would not be possible to invert for 
d'f unless the geometry were quite specific. Additionally, stability is lost since we have taken 
the Laplacian of the data. However, we will see that by taking more background solutions we 
regain ellipticity and stability. If n > 3, consider the n background solutions {ui} together. 
We in fact have that 

n 

^AdHi = (n- 2)A, 

i=l 

which is clearly elliptic, and can be inverted to find ^7 if we impose Dirichlet boundary 
conditions. Therefore, we can obtain the stable reconstruction of (^7, 6a) by inverting a 
second order system of equations. 

In dimension n = 2, the measurements Hi and H2 alone are not quite enough since 

Ad^i = -AdH2. 

We need the contribution 

H12 = Vui • VU2 

which can be obtained by the background solution = ui + U2, and yields 

AdHi2 = —2dxiX2- 

We observe that we now want to solve the 2x1 system 

-2dx,x, / \ AdHi26-f 
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One can approach this in the same manner as above and apply Aj to both sides. This yields 
the simple bi-Laplace equation 

^* ( ^dHuS^ \ 
AA(57 = A*j\ 

\ AdHi26-/ J 

which we can clearly invert stably if we impose Cauchy data 5^ = 0, = on dX. Again we 
have converted the system to one that is square with the cost of having to go to fourth order 
equations and requiring more boundary data. Note that by continuity, the leading operators 
will also be invertible for (7,0") in the vicinity of (1,0) (and indeed in the vicinity of any 
positive constants). 
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